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Problem  Studied 


This  grant  covered,  in  part,  the  cost  of  one  graduate  student  working 
on  the  NMR  imaging  of  solids  in  conjunction  with  research  sponsored  by 
the  US  ONR  under  grants  N00014-86-K-0149  and  N00014-86-C-0724  subcontract 
SK80D7820F  from  Lockheed  corporation. 

The  goal  of  the  research  was  to: 

i)  Improve  the  quality  of  NMR  images  obtained  from  simple  inorganic 
solids  so  they  are  useful  non-destructive,  non-invasive  tools  to  study 
mechanically  induced  defects  in  solids. 

ii)  Use  the  techniques  on  mechanically  damaged  materials  to  establish 
the  utility  of  the  technique  to  study  mechanically  damaged  materials. 


Summary  of  Most  Important  Results 

O  O 

The  quality  of  images  using  the  Na  signal  from  NaCl  in  a  propellant  simulant 
has  improved  as  a  function  of  time,  as  is  illustrated  in  figure  1.  The  image 
labeled  "1989"  is,  in  fact,  of  high  enough  resolution  that  the  use  of  contour 
plots  is  misleading.  The  images  are  normally  viewed  on  the  computer  screen 
however  a  simple  cost-effective  means  of  reporducing  these  images  has  not  been 
found.  The  major  improvements  to  the  equipment  are  the  increase  in  gradient 
strength  used  to  form  the  image  and  some  improvements  in  the  software  data 
processing.  The  details  of  these  improvements  has  been  previously  reported 
(see  semi-annula  progress  reports,  this  grant). 

The  contents  of  the  propellant  simulant  are  shown  in  table  I.  The  improvements 
in  our  instrumentation  were  sufficient  that  we  were  able  to  obtain  a  usable 
image  using  the  Aluminum  signal  from  the  metal  particles  in  the  simulant.  This 
image  is  shown  in  figure  2.  The  aluminum  NMR  signal  is  several  times  broader 
than  the  Na  signal  making  images  that  much  more  difficult  to  obtain.  To  the 
best  of  our  knowledge,  no  one  has  obtained  a  better  image  using  aluminum  with 
any  technique.  Attempts  at  Lockheed  corporation  to  use  sophisticated  multiple 
pulse  techniques  did  not  yield  an  improved  resolution  due  to  significant 
distortions  from  off  resonance  signals. 

The  utility  of  the  NMR  imaging  technique  to  study  damaged  materials  was 
demonstrated  and  is  contained  in  the  reprint  which  is  appendix  A  of  this  report. 
It  is  important  to  note  that  the  line  narrowing  techniques  commonly  used 
for  solids  work  cannot  be  used  to  study  quadrupole  broadened  NMR  signals 
to  study  defects.  The  line  narrowing  techniques  reduce  or  remove  the 
effects  of  the  defects  on  the  NMR  signal.  Our  work  so  far  has  concentrated 
on  the  non-destructive  measurement  of  defect  densities  in  mechanically 
damaged  NaCl  using  the  NMR  signal  from  the  Na  nucleus.  The  results  demonstrate 
the  utility  of  the  technique  to  study  mechanically  damaged  materials  using 
the  non-invasive  NMR  imaging  technique. 
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Nuclear  magnetic  resonance  imaging  as  a  tool  to  characterize  defect 
densities  in  solids 

B.  H.  Suits  and  J.  L.  Lutz 

Physics  Department,  Michigan  Technological  University.  Houghton.  Michigan  49931 

(Received  23  September  1988;  accepted  for  publication  13  January  1989) 

Nuclear  magnetic  resonance  imaging  methods  demonstrate  the  spatial  variation  of  21Na  signal 
intensity  in  NaCI  which  results  from  an  impact-induced  inhomogeneous  defect  distribution. 


Nuclear  magnetic  resonance  (NMR)  imaging  of  solids 
is  a  relatively  new,  nondestructive,  noninvasive  measure¬ 
ment  technique  which  has  been  shown  to  be  useful  in  moni¬ 
toring  inhomogeneities  of  chemical  species, 1,2  and  for  identi¬ 
fying  cracks  in  materials.1-3  The  work  presented  here 
demonstrates  the  utility  of  the  technique  to  study  the  spatial 
distributions  of  other  types  of  material  defects  as  well,  in 
particular,  dislocations  in  solids  containing  a  nucleus  with  a 
nonzero  quadrupole  moment. 

The  use  of  NMR  to  study  mechanically  induced  defects 
is  not  new.4  A  comprehensive  review  of  the  effects  on  the 
quadrupole  interaction  due  to  point  defects  and  dislocations 
in  cubic  systems  has  been  presented  by  Kanert  and  Mehr- 
ing.5  What  is  new  here  is  the  application  of  these  ideas  in  a 
spatially  resolved  mode. 

This  discussion  is  limited  to  the  use  of  the  electric  quad¬ 
rupole  interaction,  though  other  interactions  which  are 
modified  in  the  presence  of  defects  may  prove  useful  as  well. 
Thus,  the  discussion  is  limited  to  nuclei  with  spin,  />  1/2. 

The  electric  quadrupole  interaction  is  due  to  the  pres¬ 
ence  of  an  electric  field  gradient  (efg)  at  the  nucleus  under 
study  which  couples  to  the  nuclear  electric  quadrupole  mo¬ 
ment  ( eQ ) .  The  quadrupole  interaction  is  discussed  in  detail 
elsewhere.6-®  The  efg  can  be  described  using  a  symmetric, 
traceless,  second-rank  tensor  V  with  components 
Va0  (aji  =  xj'.z).  The  tensor  represented  using  one  coordi¬ 
nate  system  (x,^^)  can  be  transformed  to  another 
(xj  ,x(  jc\ )  using  the  rotation 

=  t~'Vu,T,  (1) 

where  the  elements  of  T,  Ttj,  are  the  cosines  of  the  angles 
between  x‘  and  x;.  It  is  convenient  to  discuss  the  efg  in  a 
reference  frame  where  Va0  =  0  if  a^0,  known  as  the  princi¬ 
ple  axis  system. 

For  NMR  measurements  in  a  high  magnetic  field  H,  it  is 
often  the  case  that  the  Zeeman  interaction,  described  by  the 
Hamiltonian  =  —  yft H*I,  where  y  is  the  gyromagnetic 
ratio  and  I  is  a  spin  operator,  is  large  compared  to  the  quad¬ 
rupole  interaction,  3fQ.  If  the  convention  is  used  that  the 
magnetic  field  H  is  along  the  z  axis  of  the  lab  frame  (x,y,z) 
and  the  principle  axes  of  q  define  a  crystal  reference 
frame  (.x'yj),  then  (to  first  order  in  perturbation  theory) 
one  has  a  total  interaction  described  by  the  Hamiltonian 

*  =  -  +  ( eQ /4/( 21-  I)][3/2  -  /2]  Va, 

(2) 

where,  using  Eq.  (1), 


Va  —  V2  I  cos2  8  +  Vxx  sin2  8  cos2  8 

4-  V/y.  sin2  8 sin2  <f>,  (3) 

where  6  and  <f>  are  the  Euler  angles  relating  the  primed  to  the 
unprimed  coordinates.  This  yields  magnetic  dipole  transi¬ 
tion  frequencies  com  between  the  states  characterized  by  m 
and  m+  l,m=  —  7,  —  7  +  1,...,/, 

com=yH+  [3*0/87(2/ -  l)*](2/w  +  \)VU.  (4) 

Note  that  the  transition  from  m  =  —  J  to  m  =  +  ( (present 
for  nuclei  with  half-integer  spin )  is  not  sensitive  in  first  order 
to  the  presence  of  the  quadrupole  interaction,  while  other 
transitions  can  be  strongly  affected. 

There  are  many  types  of  defects  in  materials.®-10  Asso¬ 
ciated  with  any  defect  is  a  strain  field,  the  nature  of  which 
depends  on  the  particular  type  of  defect.  The  distance  over 
which  these  strain  fields  have  an  impact  on  the  NMR  signal 
may  be  quite  large. 

For  simplicity,  the  details  of  only  one  type  of  defect,  an 
isolated  dislocation,  will  be  presented  here.  The  features 
which  emerge  for  others  are  qualitatively  similar. 

It  is  quite  difficult  to  treat  nuclei  which  are  very  close  to 
the  defect,  but  outside  this  “core”  region  elastic  theory  can 
be  used.  Estimates  of  the  size  of  the  core  region  show  that  it  is 
only  a  few  lattice  constants  in  radius.10 

In  the  region  where  the  elastic  theory  applies  the  atoms 
will  be  displaced  only  a  small  amount  relative  to  their  neigh¬ 
bors  and  hence  one  can  assume  that  the  changes  in  the  quad¬ 
rupole  interaction  will  be  linearly  related  to  the  local  strain, 
e,  or  stress,  a.  With  this  assumption  one  can  write 

K,  =  V°,  +  X  Cijklakl.  (5) 

k.t 

- 

The  coefficients  Cljkl  are  measured  quantities  and  values  for 
many  cubic  systems  are  tabulated  in  Ref.  5. 

A  dislocation  is  characterized  by  a  line  through  the  ma¬ 
terial  and  a  slip  plane.  If  the  dislocation  line  defines  the  z" 
axis  and  the  slip  plane  lies  in  the  x"  —  z '  plane,  then  the 
stress  tensor  for  a  screw  dislocation  is 

(0  0  -  sin  6\ 

0  0  cos0  )  (6) 

—  sin  8  cos  8  0  / 

and  for  the  edge  dislocation 
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where  G  is  the  shear  modulus,  b  is  the  magnitude  of  the 
Burgers  vectors,  D  —  Gb/[2ir(  1  —  v)J  where  vis  the  Pois¬ 
son  ratio  and  the  r  and  9  are  the  positions  in  the  sample 
relative  to  the  dislocation  in  the  dislocation’s  reference 
frame.  Thus,  the  efg  in  the  lab  reference  frame  can  be  calcu¬ 
lated  as 

V=T-'U-\V°  +  Ca)UT,  (8) 

where  U  transforms  the  interaction  from  the  double  prime  to 
prime  coordinates  and  Tinto  the  laboratory  frame. 

It  can  be  seen  that  both  types  of  dislocations  yield  an 
angularly  dependent  response.  In  some  cases  this  angular 
dependence  can  be  used  as  a  means  of  characterizing  the 
defects  which  are  present. 

Due  to  the  cubic  symmetry  at  the  23Na  site  in  NaCI  the 
efg  tensor  in  the  undistorted  lattice,  V°,  is  zero  and  all  the 
transitions  occur  at  the  same  frequency,  yH /fi.  As  an  explic¬ 
it  example,  the  screw  dislocation  in  NaCI  [  ( 1 10)  ( 1 10)  slip 
system]  gives 

Vzl  =  ( GbC , ,  /2rrr)f  ( 9,  T,  U) . 

The  function  f(9,T,U)  (see  Ref.  5)  describes  the  angular 
dependence  and  in  this  case  has  the  property  that  it  is  zero  if 
the  magnetic  field  is  along  or  perpendicular  to  the  disloca¬ 
tion.  Putting  in  numbers  one  has  an  additional  splitting  of 
the  /=  +  j'  to  +  i  transition  given  approximately  by 

Av=  ±  (2AUHz/r)f{9,T,U),  (9) 


of  the  spectrometer  (in  our  case  roughly  0.3  MHz).  The 
signal  from  these  regions  will  be  reduced  in  two  ways:  ( 1 ) 

the  +^to  ±  ( transitions  will  be  missing  and  (2)  the  _ i  to 

+  ]  transition  will  be  decreased  in  magnitude  due  to  the  fact 
that  it  must  be  treated  in  the  effective  spin  ( approximation.'’ 
Outside  this  region  is  a  region  labeled  “in  baseline.”  From  a 
simple  point  of  view  one  regards  the  satellite  transition  fre¬ 
quencies  from  this  region  to  be  spread  out  over  such  a  wide 
range  that  the  total  response  within  the  spectrometer  band¬ 
width  is  effectively  flat.  These  signals  are  observed  in  a  rigor¬ 
ous  sdnse  but  if  the  baseline  is  corrected  before  or  during 
image  processing,  they  are  removed.  All  that  is  observed 
then  is  the  —  ( to  ( transition. 

Outside  all  these  regions  the  shift  is  small  compared  to 
the  NMR  linewidth  (in  this  case  roughly  3  kHz)  and  hence 
the  effects  of  the  dislocation  are  obscured  and  the  full  inten¬ 
sity  is  observed. 

Of  course,  the  dividing  lines  between  the  regions  are  not 
sharp  and  there  is  some  change  in  the  central  lineshape  used 
to  create  the  image.  The  changes  in  lineshape  are  not  over¬ 
whelming  and  this  simple  picture  reflects  the  essence  of  the 
behavior. 

The  sensitivity  of  the  NMR  technique  is  quite  good.  Us¬ 
ing  the  estimated  100  nm  as  the  maximum  range  of  the  ef¬ 
fects  of  the  dislocation  then  in  a  cross  section  one  has  rotigli- 


where  r  is  the  distance  from  the  dislocation  line  in  ang¬ 
stroms.  Figure  1  illustrates  the  effects  of  the  dislocation  on 
the  NMR  signal.  Curves  of  constant  frequency  shift  of  the 
±  2 t0  ±i  satellite  transition  appear  as  lobes  about  the  dis¬ 
location. 

In  the  core  region  one  does  not  know  the  efg  but  the 
number  of  nuclei  involved  is  quite  small  and  this  region  can 
be  ignored.  Outside  the  core  region  is  a  region  where  the 
quadrupole  shift  moves  the  satellites  outside  the  bandwidth 


FIG.  1 .  A  schematic  diagram  of  the  spatial  dependence  of  the  contribution 
to  the  NMR  signal  from  the  ±  j  to  ±  J  satellite  transitions  near  a  disloca¬ 
tion.  The  figure  is  scaled  appropriately  for  the  ”Na  signal  from  NaCI  with 
the  magnetic  field  along  ( 100). 


I -  I  cm -  | 

FIG.  2.  NMR  images  showing  thespatia!  distribution  ofthe  NMR  intensity 
(a)  before  and  (b)  after  extensive  damage  to  part  ofthe  crystal  by  impact. 
The  difference  (percent)  between  (a)  and  (b)  is  shown  in  (c).  The  change 
in  signal  intensity  is  due  to  the  presence  of  defects  which  alter  the  NMR 
response.  The  sample  is  I  X0.54  cm2  and  0.2  cm  thick.  The  view  is  through 
the  thickness  of  the  sample. 
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ly  1.6x10  111  cm*  which  Is  unobserved  for  each  dislocation. 
Thus,  dislocation  densities  larger  than  about  10h  or  107/cm2 
should  be  observable  in  NaCl  using  this  technique  for  mea¬ 
surements  with  a  good  signal  to  noise  ratio. 

Figure  2  shows  images  for  a  crystal  which  has  been  da¬ 
maged  by  impact  using  a  drop  hammer  (495  g  from  9  cm 
high).  These  images  are  two-dimensional  projections 
through  the  thickness  of  the  sample  with  the  strike  along  the 
direction  of  view.  Part  of  the  sample  ( the  right-most  quarter 
in  Fig.  2 )  was  outside  the  region  of  the  strike.  The  damage  is 
clearly  evident  as  a  large  decrease  in  intensity  in  the  region 
that  was  hit.  Plastic  deformation  of  the  damaged  region  is 
roughly  5%  (there  was  no  measureable  loss  of  mass). 

Figure  3  shows  cross  sections  through  the  center  for  two 
orientations  of  the  damaged  crystal  with  respect  to  the  ap¬ 
plied  magnetic  field.  The  view  of  the  sample  is  the  same  in 
both  cases.  There  is  a  significant  difference  which  is  not  uni¬ 
form  from  one  side  of  the  sample  to  the  other.  That  is,  the 
defect  distribution  as  seen  by  the  2JNa  nuclei  does  not  have 
the  cubic  symmetry  of  the  host  lattice  and  the  degree  to 
which  the  symmetry  is  broken  varies  across  the  sample.  One 
w'ay  the  symmetry  can  be  broken  is  if  there  are  more  disloca¬ 
tions  along  one  direction  than  another  equivalent  direction. 
In  cases  where  the  defect  density  is  high  and  a  nucleus  can 
“see”  more  than  one  at  a  lime  (for  NaCl  a  local  density  of 
roughly  lO'Vcm  or  greater),  a  difference  in  the  microsco¬ 
pies  of  the  distribution  for  different  directions  may  also  give 
rise  to  this  effect.  As  a  simple  example,  if  in  one  direction  the 
dislocations  occur  in  bunches,  whereas  in  another  they  are 
spread  out  uniformly,  one  can  expect  this  phenomena.  It 
should  be  noted  that  the  original  sample,  which  was  cleaved 
from  a  larger  melt-grown  block,  also  showed  a  lack  of  cubic 
symmetry  for  the  NMR  response,  however,  in  this  case  the 
change  in  signal  intensity  upon  rotation  of  the  magnetic  field 
was  uniform  across  the  sample.  The  observed  lack  of  cubic 
symmetry  in  the  undamaged  sample  is  not  unexpected  based 
on  previous  results  for  NaCl.5 

The  NMR  data  presented  in  this  paper  was  obtained 
using  a  “homemade”  microcomputer  controlled  NMR  spec¬ 
trometer  with  a  magnetic  field  of  approximately  85  kG  (360 
MHz  for  protons,  95  MHz  for  2,Na).  The  projection  recon¬ 
struction  (convolution  backprojection )  algorithm  described 
in  Ref.  1  was  used  to  create  images  from  Fourier  trans¬ 
formed  free  induction  decays.  The  gradient  strength  used  for 
images  was  approximately  15  G/cm.  In  all  of  these  images, 
the  linear  resolution  is  limited  by  the  NMR  linewidth  to 
about  1  mm. 

The  imaging  method  used  so  far  in  these  studies  has 
been  the  bare  minimum.  With  improvements  in  the  tech¬ 
nique  one  can  hope  to  improve  the  sensitivity  and  resolution 
of  the  images  and  to  streamline  the  analysis  of  the  defects 
present. 

These  results  show  that  NMR  imaging  can  be  used  to 
measure,  and  to  some  extent  characterize,  defect  distribu- 


F1G.  3.  Cross  sections  through  images  of  the  heavily  damaged  crystal  for 
two  orthogonal  orientations  of  the  magnetic  field;  magnetic  field  (a)  per¬ 
pendicular  and  (b)  parallel  to  the  sample  face.  The  defect  distribution  does 
not  have  the  cubic  symmetry  of  the  NaCl  lattice  and  hence  destroys  the 
cubic  symmetry  of  the  NMR  response. 


tions  in  solids.  The  NMR  imaging  technique  complements 
existing  techniques  which  allow  one  to  visualize  individual 
defects.  The  advantage  of  NMR  imaging  over  these  tech¬ 
niques  is  that  one  can  see  into  the  interior  of  a  sample,  inde¬ 
pendent  of  the  material’s  opacity,  in  a  noninvasive  way.  In  an 
industrial  setting,  such  as  for  the  manufacture  of  explosives, 
NMR  imaging  has  some  obvious  advantages. 
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